Spin-dependent tunneling through a symmetric semiconductor 

barrier: the Dresselhaus effect 



L. G. Wang, Wen Yang and Kai Changl^ 

NLSM, Institute of Semiconductors, Chinese Academy of Sciences, 
P. 0. Box 912, Beijmg 100083, China 

K. S. Chan 

Department of Physics and Materials Science, 
City University of Hong Kong, Kowloon, Hong Kong, China 

Abstract 

Spin-dependent tunneling through a symmetric semiconductor barrier is studied including the 
Dresselhaus effect. The spin-dependent transmission of electron can be obtained analytically. 
By comparing with previous work(Phys. Rev. B 67. R201304 (2003) and Phys. Rev. Lett. 93. 
056601 (2004)), it is shown that the spin polarization and interface current are changed significantly 
by including the off-diagonal elements in the current operator, and can be enhanced considerably 
by the Dresselhaus effect in the contact regions. 
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Recently electron spin in semiconductors has attracted a rapidly growing interest due 
to its potential application in spintronics devices. To successfully incorporate spin into 
existing semiconductor technology, one has to overcome technical difficulties such as efficient 
spin-polarized injection, transport, control and manipulation, as well as measurement of 
spin polarization. The injection of spin-polarized electrons from ferromagnetic metals into 
semiconductors has low efficiency, less than 1%, because of a large resistivity mismatch 
Between fe™.a.„e«ca„d...o„d.c.o.„.«a.Q . Rashba proposed that this problem 
could be solved by inserting a tunneling barrier at the metal-semiconductor interface P]. 
Asymmetric nonmagnetic semiconductor barriers are also used in the construction of spin 
filterSfSi]- This effect is caused by the interface-induced Rashba spin-orbit couplingi43 and can 
be quite significant for resonant tunneling through asymmetric double-barrier structure j^. 
Very recently, a multichannel field-effect spin-barrier selector was investigated theoretically 

spin polarization and an 



interesting "tunneling spin-galvanic" effect were found in the tunneling of_e 
a single symmetric barrier utilizing the Dresselhaus effect in the barrier^ 



ectron through 
However, the 

off-diagonal elements in the current operator, the contribution from the Dresselhaus effect, 
are neglected in these works. These off-diagonal element in the currents operator could lead 
to the significant correction of the spin-dependent transmission, especially in a thin barrier 
case. 

In this paper, we investigate theoretically the spin-dependent tunneling through a single 
symmetric barrier. The barrier and contacts consist of a zinc-blende-structure semiconductor 
lacking the inversion symmetry. It is shown that the spin-dependent tunneling and the 
electric current in the plane of the interfaces are different from the previous studies 3] ■ 
This difference is significant in thin barrier case and disappears gradually with increasing 
the thickness of the barrier. It is interesting to notice that the spin polarization and the 
interface current j\\ are enhanced by including the Dresselhaus effect in the contact regions. 

We consider the transmission of an electron with initial wave vector k = (fc,;, kz) through 
a flat potential barrier of height V grown along the z [001] direction; fc,, is the in-plane wave 
vector, and kz is the wave vector along the growth direction, i.e., z axis. Then the electron 
Hamiltonian including the spin- dependent k^ Dresselhaus term is 
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H 

Hd 



2m* 



+ V{z) + Hd, 



li[(TxK{kl - kl) + ayky{kl - kl) + a-,k^{kl - /cj)] 



(1) 



where m* is the effective mass of electron, and V{z) is the height of barrier. describes 
the Dresselhaus spin-orbit couphng, a a are the Pauh matrices, and 7j(i = 1,2) describe the 
strength of the Dresselhaus effect in the contact regions and the barrier. We assume that 
the kinetic energy of the electron is substantially smaller than the barrier height V in the 



barrier 



m 



(we take V = 0.2eV, Ep = 0.02eV in our paper). The Hamiltonian in the 



barrier is simplified to 



H 



p2 Q2 

— - + V{z) + 72(fcx<7^ - h^y)^ 



2m 



The eigenvalues and eigenstates of the Hamiltonian are 

n\ki + qi) 



E 

x± 



2m* 

' ( ' 



(2) 
(3) 

n\ P (4) 

which correspond to the "+" and "-" electron states with opposite spin orientation. The 
wave functions of the electron in the left source, the barrier and right drain are 



^ffe = exp(ifc„ ■ p)2_^[Aj exp{qjz) + Bj exp{-qjz)]xj, 



= exp(ifc„ ■ p)'^tj exp{ikjz)xj, 



(5) 



here = exp(ik„ ■ p) exp {ikj z) Xj corresponds to the injected spin state of electron^, 0], 

3=± 

t±, are the transmission and reflection coefficients for the spin states x±; respectively. 
p = {x, y) is the coordinate in the barrier plane, and ip is the polar angle of the wave vector 
k in the xy plane. The wave vectors q±{k±) inside and outside the barrier are given by 



9± 



\ 1 ± ¥72A:„ 



(6) 
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where A = a/1 + tan^ 6 + tan^ 6 sin^ 2y9(tan^ 9/ A — 2) and mKm'^) is the effective masses 
outside(inside) the barrier. The boundary conditions inchiding the Dresselhaus effect are 





z=0 


= ^b\z=0, 




z=a 




JL^L 


z=0 


= jb'^b\z=0, 


Jb^b 


z=a 





(7) 

where = L, R, b) are the current operators in the left side, the right side of the barrier, 
and inside the barrier. 



Jl,r — j: 



' dz ' m'y dz 



4 -272fc„e^^ \ d 

Note that the off-diagonal elements of the current operators were neglected in Ref.j3]- For 
the real case ml'yik„/h^ <^ 1 and ml'y ikf^cos{2(f )/q±/h^ <^ 1, we obtain the analytical 
expression of the spin-dependent transmission coefficient of the electron 

f, = ^ ^ ^^"^ (Q) 

where 6± = 1 ± 2mljik„/ h'^ , and 6'^ = 1 ± 2m*2'^2ku/fi^- 

It is convenient to introduce the spin polarization P determined by the difference between 
the transmission of the spin states x+ and X-, 



P = '-hHtt^^ (10) 
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The interface current due to spin-polarized electron transport through the tunneling 
structure can be obtained through the spin density matrix ^ 

•^'ll = ^ 5Z TpTr[TpiT+v^]v,,, (11) 

fc„,fcz>0 

where Tp is the momentum relaxation time, pi is the electron density matrix on the left 
side of the structure, and T is the spin matrix of the tunneling transmission that links the 
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incident spin wave function to the transmitted spin wave function "^r, "^r — T^l. The 
spin matrix of the electron transmission through the structure is given by 



r^Y^^sXsxt- (12) 



s=± 

In the case of small degree of spin polarization, the density matrix has the form 

P. = /o/-^^^(n.-a), (13) 

where /o is the equilibrium distribution function of non-polarized carriers, Ps is the degree of 
the polarization, and < 1/s > is the average of the reciprocal kinetic energy of the carriers. 
Us is the unit vector directed along the spin orientation, and the orientations of spin s± in 
the states "+" and "-" depend on the in-plane wave vector of the electron and are given by 
s± = (=F cos ± sin 0). Taking into account the spin matrix and the density matrix, the 
interface current is 



n--Co E - \U']v.v^,. (14) 



dfo 

k„,kz>0 

where Co = er^ ^^^^ ■ It is interesting to notice that the direction of the interface current 
depends on the spin polarization of the injected electron pg. 

It is well know that normally the transmission through the barrier reaches maximum for 
carriers propagating along the normal to the barrier in the absence of spin-orbit interaction. 
But the spin-orbit coupling changes this rule as shown in Fig. l(see Eq. (9)). The tunneling 
transmission for the spin-polarized electron with the finite in-plane wave vector k\\ is larger 
than that for the electron with the opposite in-plane wave vector, —fey. This asymmetry 
results in the in- plane flow of the transmitted electron near the barrier, i.e., an interface 
electric current. 

For an electron with wave vector fe(see in Eq. (5)) is injected into the barrier, the 
spin-dependent transmission is determined by the width of the barrier and the Dresselhaus 
spin-orbit coupling strength(see Eq. (8)). The spin-dependent transmission and the spin 
polarization P are plotted as a function of the barrier width q^a in Fig. 2. Here we consider 
the Dresselhaus effect only in the barrier material. The material parameter relevant to GaSb 
are used in our calculation 71 = 0, 72 = 187 and the effective masses are m* — rn^ — 0.041mo. 
It is shown that the spin-dependent transmission decreases rapidly with increasing barrier 
width goo, while the spin polarization increases gradually. From this figure, we can see 
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that the results of Ref . [7[ agree with ours only in the thick barrier case, but there is a big 
difference in the thin barrier case. From the numerical results, we can find that the off- 
diagonal elements due to the Dresselhaus effect, in the current operator play an important 
role in the spin-dependent transport, especially in the thin barrier case. 

Figure 3 describes the dependence of the magnitude of the interface current j\\ on the 
width of the barrier in unit of Cq for incident electrons, which form spin-polarized degenerate 
gas as in GaAs(Fig. 3(a)) and in GaSb(Fig. 3(b)), respectively. From this figure, we can 
see that the two results are very different. Our numerical results show that the maximum 

with increasing the barrier width. The numerical accuracy spin-dependent tunneling and 
interface current could be improved significantly by including the effect of the off-diagonal 
elements of the current operator due to the Dresselhaus effect. The difference between the 
present results and Ref.jsl will increase with increasing 7712^1,72/^^ comparing Fig. 3(a) and 
3(b). 

The spin polarization P is plotted as a function of the barrier width a for different 
strengths of the Dresselhaus effect in Fig. 4. From this figure, we see that the spin polariza- 
tion depends on the Dresselhaus effect 7j(? = 1,2) and the barrier width a. It is interesting 
to notice that the spin polarization is enhanced by including the Dresselhaus effect 71 in 
the contact regions, and increases with increasing the barrier width a. On the other hand, 
the spin polarization decreases with diminishing the strength of the Dresselhaus effect 72 
in the barrier for a fixed 71. The spin polarization saturates gradually with increasing bar- 
rier width a for 72 = 0, and reaches maximum when 71 = 72 = 187. These features can 
be understood from Fig. 4(b) and 4(c). The total transmission decreases with increasing 
barrier width and almost the same for different strengths of the Dresselhaus effect (see Fig. 
4(c)), but the differences between the spin states x+ and X- reach the maxima for different 
strengths 71, 72 (see Fig. 4(b)). Therefore the spin polarizations exhibit different behavior 
for different strengths of Dresselhaus effect in the contact and barrier regions. The spin 
polarization saturates gradually with increasing barrier width for 71 = 187, 72 = since the 
spin polarization is mainly caused by the Dresselhaus effect 71 at both sides of the barrier, 
and consequently almost independent of the barrier width. The spin polarization approaches 
the maximum when 71 = 72. 

Fig. 5 shows the interface current jii as a function of the barrier width a for different 
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strengths of the Dresselhaus effect. Prom this figure, we can see that the interface current 
j|l is enhanced by the strength of the Dresselhaus effect 71 in the contact regions, and show 
maxima at certain values of the barrier width a. It is interesting to notice that the interface 
current jy reaches maximum when 71 = 72, and decreases as the strength of the Dresselhaus 
effect either in the barrier or in the contact regions decreases. Considering the kinetic energy 
of the electron is substantially smaller than the barrier height V and e~*±" <^ 1, we can 
obtain 

Jll oc - oc ((5;(5+e-'«+" - S'J^e-'"'-''). (15) 

It is obvious that the interface current j\\ increases with increasing 72 if we fix 71. Similarly, 
the interface current j\\ decreases with decreasing 71 when we fix 72. Therefore the interface 
current j\\ reaches a maximum when 71 = 72. 

In conclusion, we investigated the spin-dependent tunneling through a symmetric semi- 
conductor barrier consisting of zinc-blende semiconductor material. It is interesting to notice 
that the spin polarization and the interface current jy are enhanced significantly by includ- 
ing the Dresselhaus effect in the contact regions, and they all reach a maximum when the 
strength of the Dresselhaus effect in the barrier is equal to that in the contact regions, i.e., 
7i = 72- 
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FIG. 1: The spin-dependent transmission of electron as a function of the angle 6 between the 
direction of k and the z axis. 



FIG. 2: The spin-dependent transmission of electron as a function of barrier width qoa for GaSb 
material(7i = 0, 72 = 187, V = 0.2eV). The inset shows the spin polarization P as a function of 
barrier width qoa for the same condition. 



FIG. 3: The interface current j\\ as a function of barrier width a in unit of Co for (a) GaAs(7i = 
0, 72 = 24,V = 0.2eV) and (b) GaSb(7i = 0, 72 = 187, V = 0.2eV). The insets show the difference 
in transmission between the spin states x+ a^^d 
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FIG. 4: (a) The spin polarization of electron as a function of barrier width a including the Dres- 
selhaus effects inside and outside the barrier, (b) The difference between the transmission of the 
spin states x+ S'ld %_ . (c) The total of the transmission of the spin states x+ s^nd %_ . 

FIG. 5: The interface current as a function of barrier width a in unit of Co including the 
Dresselhaus effects inside and outside the barrier. 
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